The Superfluid Density in the D-density-wave Scenario 
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Recently Chakravarty, Laughhn, Morr and Nayak made an interesting proposal that the cuprate 
superconductors possess a hidden "d-density wave" order. We study the implication of this proposal 
for the superfluid density ps. We find that it predicts a temperature gradient \dps / dT\T=(i that is 
strongly doping dependent near the critical doping at which the superconducting gap vanishes. 



PACS numbers: 74.25.Jb, 79.60.-i, 71.27.+a 

Recently Tallon and Loram critically examined the ex- 
isting experimental data on specific heat, photoemission, 
magnetic susceptibility and optical conductivity of the 
cuprate superconductors They argued that these 

data are consistent with the existence of two competing 
energy gaps, a pseudogap and the superconducting gap, 
both of which disperse according to d-wave symmetry. 
The pseudogap exists in the normal state of the under- 
doped cuprates while the superconducting gap opens up 
in the BCS fashion at temperature Tc- Tallon and Lo- 
ram also argue that the data indicates the closing of the 
pseudogap at a doping level x ~ 19%. 

This observation motivated Chakravarty, Laughlin, 
Morr and Nayak (CLMN) to propose the "d-density- 
wave" (DDW) scenario ||^. According to this proposal 
a staggered long-range order in the orbital magnetic mo- 
ment, as firstly suggested by Hsu, Marston and Affleck in 
another context, exists at low temperatures when the 
doping is less than a critical value Xc ~ 0.19. Such order 
gives rise to a non-superconducting gap (the pseudogap) 
with the same (cosfc^; — cosfcj,) dispersion as the super- 
conducting (SC) gap. Moreover, the DDW order com- 
petes with superconductivity and causes Tc{x) to trace 
out the familiar "superconducting dome" in the doping- 
temperature phase diagram of the cuprates. So far there 
is no direct experimental evidence for the DDW order. 
Neither is there evidence for the expected Ising-like phase 
transition into such a state. CLMN argue that disorder 
gets in the way of a sharp phase transition and turns it 
into a cross-over. 

In addition to offering a possible new phase for high- 
Tc cuprates, there are a number of attractive features in 
CLMN's proposal. Since the DDW metal (i.e. the normal 
state in the presence of DDW order) is a doped band- 
insulator one expects the Drude weight in the normal 
state , and the superfluid density in the superconduct- 
ing state to be both proportional to the doping density 
as indeed found in the experiments. The closing of 
the DDW gap at x = 0.19 can easily explain the observed 
kink in the jump of the T-linear specific heat coefficient 
In addition the small hole pockets centered around 
the nodes of the DDW gap could be the progenitor of the 



Fermi arcs observed in angle- resolved photoemission [Q. 
From a technical point of view a doped band insulator 
has more resemblance to a doped Mott insulator than a 
large-Fermi-surface metal, hence is a better starting point 
for description of the underdoped cuprates. Finally the 
CLMN proposal is crisp and in principle falsifiable. For 
all the above reasons we feel that it is worthwhile to fur- 
ther check the prediction of this proposal against existing 
experiments. 

The theory presented in Ref. is mean-field in na- 
ture. For such a theory, quasiparticles are sharply de- 
fined. Since a much better case can be made for the exis- 
tence of quasiparticles in the superconducting states [|j , 
we confine our calculations to low temperatures where 
SC order parameter exist. 

Among various superconducting properties we focus on 
the superfluid density ps because of its rather unconven- 
tional doping and temperature dependence, which has 
been the focus of many theories |^,^. Experimentally 
it is established that, for a fairly wide range of doping, 
dps/dT is nearly doping independent at low temperatures 
]lO| , pl| . In contrast the extrapolated zero temperature 
superfluid density changes significantly with doping . 
The goal of this paper is to work out the DDW theory's 
prediction for pa (T, x) (in the absence of disorder) . 

Even if the quasiparticles are well-defined in the under- 
doped superconducting states, the Mott constraint can 
substantially modify the results of the free, mean-field 
theory. Therefore along the mean-field prediction we also 
present the results of the "projected DDW model" where 
the electron occupation constraint is taken into account. 
More specifically, the no-double-occupancy constraint is 
implemented by introducing the slave bosons (holons) 
plus the gauge fields which couple to both bosons and 
fermions (spinous). The strict occupation constraint is 
refiected in the absence of the Maxwell term for the gauge 
fields, i.e. the coupling constant is infinity. In a recent pa- 
per, one of us looked into such a gauge theory where the 
underlying mean-field vacuum is the d-wave RVB state 
of Kotliar and Liu ||l2|. It was shown that the gauge 
field can be integrated out exactly in the continuum ap- 
proximation of the lattice theory |l3,H . This continuum 
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theory describes the (correlated) density and current fluc- 
tuations of holons and spinous above the length scale of 
the inter-holon distance A^. The physics below such a 
length scale is entirely summarized by a few parameters 
in the effective action. In the absence of more accurate 
estimates for these parameters, the author of Ref. 
took the mean-field prediction of them. 

In the following treatment of the projected DDW 
model we shall follow the same path taken in Ref. ||l^] 
while replacing the Kotliar-Liu RVB mean-field vacuum 
by the CLMN mean-field vacuum in the spinon sector. 
We expect the program carried out in Ref. to work 
well in the presence of DDW order, because the DDW 
metal-a doped band insulator-has very little spinon den- 
sity and current fluctuations below the length scale Xh- 
The fluctuatioins above this length scale is already cap- 
tured by the analysis of Ref. . 

Our results are as follows. For xi < x < Xu and at 
low temperatures, the mean-field DDW theory predicts 
a superfluid density that behaves as 



Ps{T,x) = Pddw{Q,x) - addw{x)T. 



(1) 



Results for pddw{^Tx) are shown in the main panel, Fig. 
1(a). Note that the zero-temperature superfluid density 
is non-zero at x = xi,Xu- An infinitesimal temperature 
will, however, destroy superfluidity because the pairing 
gap vanishes there. The temperature gradient of ps is 
given by addw{x) ~ t/A{x) where t is the hopping inte- 
gral and A(a;) is the maximum rf-wave superconducting 
gap. It diverges at two points, x = xi, and x — Xu 
because A(a;) vanishes there. In general the doping de- 
pendence of A (a;) is reflected in addw{x) as shown in the 
insert of Fig. 1(a). 

The projected DDW model predicts 



Ps{T,x) = Zj{x)pddw{0,x) - Zj{xYaddw{x)T, 



(2) 



where Zj (x) is the current renormalization factor to be ex- 
plained later. The doping dependence of Zj{x)pddw{^,x) 
(main panel) and Zj{x)'^adds{x) (insert) are shown in 
Fig. 1(b). Unlike the mean-field result, dps/dT has a 
sharp variation near Xc which originates from the kink 
in the mean-field superfluid density pddw{x, 0). 

Let us compare the above theoretical results with exist- 
ing data on dps (T, x) /dT. In the oxygen-depleted YBCO 
thin films reported in Ref. |lll|| , the transition tempera- 
ture Tc ranges from 90K at the optimal doping to 38K on 
the underdoped side. Meanwhile, \dps/dT\ decreases by 
about 15% from its maximal Yahie at optimal doping. In 
contrast, the mean-field DDW theory predicts an increase 
of \dps/dT\ by 140%, provided scales with Ao(x). Af- 
ter projection, the dependence of \dps/ dT\ on x is about 
100% . It should be pointed out that in the DDW theory 
\dps/dT\ is not universal, and is subject to quantitative 
change if the doping and temperature dependences of the 
pairing amplitude are different from what we assumed in 



this paper. It is not out of question that the experimen- 
tal data happens to span the plateau regime in Fig. 1(b) 
(insert). In this case, the present theory predicts a rapid 
increase of \dps/dT\ when the doping level is decreased 
even further. In the following we report the details of the 
calculations. 
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FIG. 1. (a) The doping dependence of Ps,ddii,(0, a;) (main 
pannel) and addwix) (inset), (b) The doping dependence 

of ps{Q,x) = Zj{x)pddw{Q,x) and a{x) = Zj{x)^addm{x) us- 



ing the same notation as in (a) 
Ao = 0.2t, Do = t, and h = 2t. 



The parameters used are 



The mean-field DDW theory Following CLMN, we 
adopt the following mean-field Hamiltonian 

Hddw = "^{Xk - P)cl^Cka + '^{iDkcl^Q^Cka + h.c.) 
ka ka 

- ^ Afc(cfe|C_fe| +h.c.). (3) 

k 

Here Q — (tt, tt), and Xk = '-2t{coskx + cos ky) is 
the nearest- neighbor tight-binding dispersion relation, 
Dk = D{x) (cos kx — cos ky ) is the momentum space DDW 
order parameter, and A^ = A{x){cos k^ — cosky) is the 
d-wave superconducting gap function. In the above and 
the rest of the paper we shall assume that the lattice 
constant is unity. 

Due to the breaking of translation symmetry by the 
DDW order, the Brillouin zone is half of its original 
size. For A{x) = there are two bands given by ek± = 
^^X\ -I- D\ — p In the presence of pairing, the 

quasi-particle dispersion becomes Ek± = \J ^k± + l^fcl^- 
In order to study the doping dependence of the super- 
fluid density it is necessary to specify the dependence of 
D{x) and A(a;) on x. Following CLMN, we use the so- 
lution of dE{D,A)/dD = and dE{D,A)/dA = to 
parametrize D{x) and A(x) where 



E{D, A)= aoix - Xd)D^ + aA{x - Xu)A'^ + boD^ 



-bAA-^ + wD^A^ 



(4) 



Here a^i, a a, bjj, &a, and w are (positive) material- 
dependent constants, while Xd (xu) is the doping level 
below which DDW (SC) order parameter become nonzero 
in the absence of the other. In the presence of pairing the 
onset of DDW happens at Xc < Xd while the SC order 
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parameter begins its decline precisely at Xc and vainshes ^i^) — {x — xi)/{xc — xi)\ xi < x < Xc 

at XI. The functional form for the order parameters are ^(^^ ^ Aq " 2;)/(a;„ - re,); x, < x < x„. (5) 

given by 

Here Dq (Aq) is the overall scale of D (A). In the foUow- 



D{x) = Do\/ {xc — x)/xc ing we set xi — 0.05, Xc — 0.19, Xd = 0.2 and Xu — 0.3. 

The superfluid density meausures the free energy increase caused by the phase twist of the superconducting order 
parameter |l^] . After some lengthy but straightforward algebra we obtain 

ab 1 A ^fei' 1, P^^kv \ ^ ^ P p, A 1,2 P^^ku 

Ps ( 1 " -rr- tanh — — 1 dadbSkv - -rj^ OaekvObSkv'&^ca — — 



12 ^ y Ek,. 2 ; ^ " 2L2 

k.i/—± fc,!^— lb 



^ ^ t^lAfcp tanh(/3gfc./2) {XkdgDk- DkdaXk){XkdbDk- DkdbXk) 
'L\^^El^^El_ Ek. Xl + Dl ' 



where 9a denotes the momentum derivative in the a- 
direction, L is the linear dimension of the system. 
By symmetry, the tensor is diagonal and direction- 
independent. The last term is non-zero only if the DDW 
and superconducting order coexist. It is a consequence 
of the time-reveral symmetry breaking. A similar ex- 
pression, without the last term, has been derived in Ref. 
[ psf . For the entire doping range considered, the third 
term has a negligible contribution to the superfluid den- 
sity compared to others. For self-consistency, one can 
check that the superfluid density vanishes in the absence 
of pairing due to the fact that 1) the third term van- 
ishes identically when A(a;) — 0. 2) The first two terms 
combine to give Y^ku da{[l-tan}i{f3ekv /2)]dbekv} = 0. At 
temperatures much smaller than the maximum supercon- 
ducting gap, the suppression of superfluid density comes 
from the thermal excitation of nodal quasiparticles. The 



second term in Eq. gives a T-linear suppression, while 
all the other terms contribute higher order temperature 
corrections. As a result, the superfluid density takes the 
form of Eq. (|l|) . 

We find addw{x) — (4 In 2)t/7rA(a;), independent of the 
DDW order parameter. As to Pddwi^^x) in Eq. (|l|), we 
were only able to compute it numerically, and the result 
is shown in Fig. 1(a). 

The general trend of the doping dependence of 
Pddw{Q,x) shown in Fig. 1(a) is consistent with experi- 
mental data. However the same cannot be said about 
ctddw{x). Due to the competition between DDW order 
and superconductivity in the DDW theory, A(x) is sup- 
pressed by the emergence of the DDW order, forming 
a dome. Thus A(a;) — > which implies 

that addw{x) diverg Such x-dependent 

dps/dT has not been observed experimentally. 



The projected DDW model The lattice action for the projected DDW model is given by 
L — Lb + Lddw — i aoi 

i 

Lb = Yy^i^^o+iaio-iA^- pb)bi] ~ tbYy^'^'''''^''^^^ 
Lddw = Y.^f^a{^o+ia^o-pf)ha] " ^ [(t-f i )e"- ^ /j. -f h.c] - ^[Aye'*-e<,,.^/j,. -I- h.c] (7) 



In the above bi and /io- are the holon and spinon fields 
respectively. With (pij set to zero, Lddw is the real 
space equivalent of Eq. (^). In Eq. (0) is the gauge 
field that enforces the constraint fojfei -I- fj^fia — 1- The 
DDW mean-field theory describes a doped band insula- 
tor, consequently spinon density fluctuations only occurs 
above the length scale A^j. Obviously the same is true 
for the holon density fluctuation. As a result it should 
be adequate to project out the spinon and holon density 
fluctuations above the length scale A^. 



The effective action above such cutoff length is 

H ^b ^ ^fp ^ ^Dirac ^ 

Cb = ^\ct>l{V + i^-iA)ct>b\'' + ^5pl 

^sp = ^Wsp{y + 2za)<^,pp + ^(<^:^ao0.p + 2za^f 
Cj = iSpb{(l)ldo(l)b + iao - iAo) + J^''{(t)*spd^(psp + 2ia^) 
+ p{4'bdo<Pb - ^(j)lpdo(j)sp) - ipAo - zjo ■ A. (8) 
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Here p is the doping density, and jo is the trans- 
verse ground-state current produced by the DDW or- 



der, 



and 



are the U(l) phase factors associ- 



ated with the holon field and the spinon pair-field re- 
spectively. In addition Coirac is the Dirac action for 
the spinon quasiparticles near the d-wave gap nodes, 

J'^P = ^{J2n'^la'rz'^na,iVFfpla'4'l<y,'i'VFlljlcrtp2a) is their 

3-current (r^ is the third component of the Pauli ma- 
trices, and ijjna- is the spinon Nambu spinor associated 
with the nth d-wave gap node). Kb ~ hx is the holon 
superfluid density, AKsp is the spinon zero-temperature 
superfluid density given by Eq. (m). Ub and Usp also de- 
pend on the parameters in Eq. (H), however their values 
are not important for the following discussion. 

Given Eq. (||) the gauge field can be integrated out 
straightforwardly to yield the final effective action of a 
correlated DDW superconductor 



C = 



K 



1 

2U 



- 2iA) ^ Zpp'^P{<p*dQ(l) - 2iAo) 



(9) 



In the above (j> = 4>*sp4>'bi 



K = KspKb/{Kb + AKsp), U = Usp + Wb 
Zj=Kb/{Kb + AKsp), Zp = Usp/{Usp + Wb), (10) 



and 



Di' 



2|J'?P|2 
Kb + AK, 



2{piPf 



sp 



u,- 



AUs- 



(11) 



Combining the above result with Eq. (|i|) we obtain 
Eq. (||) as the superfluid density prediction of the pro- 
jected DDW theory where 



Zj{x) = tbx/[tbX + pddwiO,x)]. 



(12) 



In Fig 1(b) we plot the x-dependence of the zero- 
temperature superfluid density in the main pannel and 
dps/dT in the insert. For x > xi the zero-temperature 
superfluid density varies with x in roughly linear fashion. 
As in the mean-field DDW theory, dps/dT has a strong 
a;-dependence. 

In conclusion, we showed that in the DDW scenario the 
zero-temperature superfluid density increases monotoni- 
cally with doping. Its temperature gradient \dps/dT\, 
being proportional to 1/A(x), also depends sensitively 
on doping. Such sensitivity to doping is expected in any 
mean-field theory where another order parameter com- 
petes with d-wave superconductivity. It is possible, how- 
ever, that this result is due to the particular approxima- 
tion scheme employed in this paper, and that another 



scheme may well render the temperature gradient insen- 
sitive to doping. 
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